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The Sudbury Neutrino Observatory (SNO) is capable of measuring simultaneously 
the flux of electron-type neutrinos and the total flux of all active flavours of neu- 
trinos originating from the Sun. A model-independent test of neutrino flavour 
transformation was performed by comparing these two measurements. Assuming 
an undistorted neutrino energy spectrum, this transformation has been definitively 
demonstrated in the pure D2O phase of the SNO experiment. In the second phase 
with dissolved NaCl in the D2O, the total active solar neutrino flux was measured 
without any assumption on the energy dependence of flavour transformation. In 
this talk, results from these measurements, their physics implications and the cur- 
rent status of the SNO experiment are presented. 



1. Introduction 

For over three decades, solar neutrino experiments^'^''^''*'^'^ have been ob- 
serving fewer neutrinos than what are predicted by the detailed models 
(e.g. *") of the Sun. This deficit is known as the Solar Neutrino Problem. 
These experiments probe different parts of the solar neutrino energy spec- 
trum, and show an energy dependence in the measured solar neutrino flux. 
These observations can be explained if the solar models are incomplete or 
neutrinos undergo flavour transformation while in transit to the Earth. 

The Sudbury Neutrino Observatory^, using 1000 tonnes of 99.92% iso- 
topically pure D2O as the neutrino target, was constructed to resolve this 
puzzle. The D2O target is contained in a 12-m diameter acrylic sphere, 
which is surrounded by 7000 tonnes of ultra-pure H2O. This volume of 
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H2O shields the D2O from high-energy 7 rays and neutrons originating 
from the cavity wall. A 17.8-m diameter stainless steel structure supports 
9456 20-cm diameter inward- facing photomultiplier tubes (PMTs). A non- 
imaging light concentrator is mounted on each PMT, extending the total 
photocathode coverage to 55%. An additional 91 PMTs are mounted facing 
outward on the support structure to serve as a cosmic-ray veto. 
The Sudbury Neutrino Observatory detects solar neutrinos by 



The charged-current (CC) reaction on the deuteron is sensitive exclusively 
to fe, and the neutral-current (NC) reaction has equal sensitivity to all 
active neutrino flavours (i/^;; x = e, 11, t). Elastic scattering (ES) on electron 
is also sensitive to all active flavours, but with reduced sensitivity to i/^ and 
Vt • By counting the free neutron in the final state of the NC reaction, the 
total active neutrino flux can be inferred for neutrinos with energy above 
the 2.2-MeV kinematic threshold. 

In the first phase of the SNO experiment, the free neutron from the 
NC interaction was detected by observing the 6.25-MeV 7 ray following 
its capture on a deuteron. To enhance the neutron detection efficiency, 2 
tonnes of salt was added in the second phase of the experiment. The free 
neutron was captured by ^^Cl and a 7-ray cascade with a total energy of 
8.6 MeV is emitted. In the final phase of the experiment, "^He proportional 
counters, the "Neutral-Current Detectors" (NCD)^, are deployed on a 1-m 
grid in the D2O. Free neutrons are captured by "^He. In this talk, solar 
neutrino results from the first two phases of the SNO experiment^^'^^'^^'^^ 
are presented, and their physics implications are discussed. 

2. Previous Results from the Pure D2O Phase 

Data from the first phase of the experiment, in which a pure D2O target 
was used, correspond to a livetime of 306.4 days. Over 450 million triggers 
were recorded during this period. 

To test the null hypothesis of no neutrino fiavour transformation, the 
extended maximum likelihood method was used to extract the CC, ES and 
neutron (i.e. NC+background) contributions in the data set. Backgrounds 
were constrained to the measured values. Data distributions in Tog-, the 
volume- weighted radial variable {R/Rav)^ = P and cos^© were simulta- 
neously fitted to the probability density functions (PDFs) generated from 



Ve + d P + P + e~ 
Vx+d p + n + V, 
Vx + e~ ^ fx + e~ 



■x 
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simulations. Rav = 600 cm is the radius of the aery he vessel, and 9q is 
the angle between the reconstructed direction of the event and the instan- 
taneous direction from the Sun to the Earth. Assuming an undistorted 
energy spectrum'^^, the fluxes of the i^e {(f>e) and non-^e {(t^^r) components 
were found to be (in units of 10^ cm^^ s^^): 

<j>^r = 3.4ltHt(«tat.)!'^;^« (syst.) 

Combining the statistical and systematic uncertainties in quadrature, (p^r 
is 5.3a above zero, and provides strong evidence for flavour transformation. 

Rcinioving the energy constraint, the signal decomposition was repeated 
using only the cos 6'© and p. The total flux of active *B neutrinos in the 
unconstrained flt was found to be (in units of 10^ cm^^ s^^) 

0„c = 6.42li j^(stat.)+°-j^ (syst.) 
3. The Salt Phase 

Neutrino flavour transformation can be an energy dependent process (e.g. 

the MSW mechanism^^). The energy constrained analysis in pure D2O 
made the explicit assumption of an undistorted *B neutrino spectrum. The 
unconstrained results suffered from the strong anti-correlation between the 
CC and NC channels. The addition of 2 tonnes of NaCl to the D2O target 
was aimed at improving the separation between CC and NC without in- 
voking the energy constraint. Neutron capture on ^^Cl typically produces 
multiple 7 rays while the CC and ES reactions produce single electrons. 
The greater isotropy of the Cherenkov light from neutron capture events 
relative to CC and ES events allows good statistical separation of the event 
types. 

Data from the salt phase that are presented here were recorded between 
July 26, 2001 and October 10, 2002, totaling 254.2 live days. The number of 
triggers recorded was over 435 million. To reduce the effects of backgrounds 
in the analysis, an effective electron kinetic energy threshold T^g > 5.5 MeV 
and a fiducial volume with radius _Rfit < 550 cm were chosen. 

Detector diagnostic triggers and instrumental background events were 
first removed in the offline data processing. The latter reduction was ac- 
complished by a set of algorithms that remove events that do not have the 
characteristics of Cherenkov light emission. For the events that passed the 
first offline reduction, the calibrated times and positions of the hit photo- 
multiplier tubes (PMTs) were used to reconstruct the vertex position and 
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the direction of the particle. The energy estimator then assigned an effective 
kinetic energy TeS to each event based on these reconstructed parameters 
and the number of hit PMTs. The energy estimator used cahbration re- 
sults such as the time variation of light attenuation at various wavelengths 
as inputs, and was normalized to the detector response to a ^^N source^^. 
^^N data taken throughout the running period verified the gain drift (~2% 
year"-*^) predicted by simulations based on the optical measurements. The 
energy scale uncertainty is 1.1%. After all the data reduction cuts and the 
selections of energy and fiducial volume, the data set was reduced to 3055 
candidate events. 

The neutron detection efficiency was calibrated with a ^^^Cf fission 
source. The detection efficiency for NC reactions in the heavy water was 
0.399 ± 0.010 (calibration) ± 0.009 (fiducial volume) for the chosen energy 
threshold and fiducial volume, an increase of approximately a factor of three 
from the pure D2O phase. This is shown in Fig. 1(a). 

Event isotropy was characterized by Pi , the mean of the Legendre poly- 
nomial Pi of the cosine of the angle between PMT hits 6ij : 



The combination /3i + 4/^4 = /3i4 was selected as the isotropy measure to 
separate NC and CC events. Systematic uncertainty on /3i4 distributions 
for simulated signal events was evaluated by comparing ^^N calibration 
data to simulations for events throughout the fiducial volume and running 
period. The uncertainty on the mean value of /3i4 is 0.87%. Comparisons of 
(3i4 distributions from '^^N events and neutron events from ^^^Cf to Monte 
Carlo calculations are shown in Fig. 1(b). The Monte Carlo calculations 
of /3i4 have been verified with a 19.8-MeV 7-ray source high-energy 
electron events (dominated by CC and ES interactions) from the pure D2O 
phase, neutron events following muons, and low-energy calibration sources. 

Radioactive decays of the daughters in the natural ^'^^Th and ^ss^j 
chains were the dominant backgrounds in the neutrino signal window. The 
levels of ^"^^Th and ^38^ daughters in the detector were determined by ra- 
diochemical assays of the detector target^^ and by analyzing the Cherenkov 
signals from their decays in the neutrino data set. The average rate of back- 
ground neutron production due to photodisintegration of the deuteron from 
radioactivity in the D2O is 0.72j^Q'23 neutrons per day. ^^Na was produced 
in the water circulation system by neutrons emitted from the rock wall in 
the underground laboratory. It could also be produced in-situ by calibra- 
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Figure 1. (a) Neutron capture efficiency versus source radial position for the pure D2O 
phase (capture on D) and salt phase (capture on CI or D) deduced from a ^^^Cf source, 
with fits to an analytic function (salt) and to a neutron diffusion model (D2O). (b) 
Event isotropy from data and Monte Carlo calculations of a -^^^Cf source and an ^''N 
7-ray source. 



tion sources (e.g. the ^^^Cf source). In order to understand its production 
and distribution, a controlled neutron activation was carried out. The am- 
plitude of this background was found to be small. 

Neutrons and 7 rays produced at the acrylic vessel and in the H2O 
can propagate into the fiducial volume. Radon progeny deposited on the 
surfaces of the vessel during construction can initiate (a,n) reactions on 
^^C, ^^O, and ^^O, and external 7 rays can photodisintegrate the deuteron 
giving rise to a free neutron. The radial distribution of these external- 
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Table 1. Background events. The internal neutron and 
7-ray backgrounds are constrained in the analysis. The ex- 
ternal-source neutrons are reported from the fit. 



Source 


Events 


Dcuteron photodisintegration 


70 i+M.U 
'•J-^-23.6 


■^H(a, a)pn 


2.8 ±0.7 


"'l»0(a,n) 


1.4 ±0.9 


Fission, atmospheric u (NC + 




sub-Chcrenkov threshold CC) 


23.0 ±7.2 


Terrestrial and reactor P's 


2.3 ±0.8 


Neutrons from rock 


< 1 


^^Na activation 


8.4 ±2.3 


n from CNO u's 


0.3 ±0.3 


Total internal neutron background 


111 

lll.J_24 g 


Internal 7 (fission, atmospheric 


5.2 ± 1.3 


l^N decays 


< 2.5 (68% CL) 


External-source neutrons (from fit) 




Cherenkov events from /3 — 7 decays 


< 14.7 (68% CL) 


"AV events" 


< 5.4 (68% CL) 



source neutrons was included in the analysis to extract this component. A 

class of backgrotmd events were reconstructed near the acrylic vessel and 
were characterized by a nearly isotropic light distribution ("AV events"). 
A limit of 5.4 events (68% CL) was set for this background. Table 1 is a 
summary of the backgrounds. 

A blind analysis procedure was used to minimize the possibility of in- 
troducing biases. The data set used during the development of the analysis 
procedures and the definition of parameters excluded an unknown fraction 
(<30%) of the final data set, included an unknown admixture of muon- 
foUowing neutron events, and included an unknown NC cross-section scaling 
factor. After fixing all analysis procedures and parameters, the blindness 
constraints were removed. The analysis was then performed on the "open" 
data set, statistically separating events into CC, NC, ES, and external- 
source neutrons using an extended maximum likelihood analysis based on 
the distributions of isotropy, cos^©, and p. The fits are shown, with sta- 
tistical uncertainties only, in Fig. 2. This analysis differs from the analysis 
reported in Sec. 2 that no assumption on the energy dependence of the 
neutrino flavour transformation was made. 

The extended maximum likeHhood analysis yielded 1339.6ttli CC, 
170.3^20;? ES, 13U.2tll[l NC, and Si.btlH external-source neutron 
events. The systematic uncertainties on derived fluxes are shown in Ta- 
ble 2. The fitted numbers of events give the equivalent fiuxes (in units 
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of 10^ cm-2s-i): 



IIP = 1.59i°:EJ?(stat);^;°^(syst) 

iSNO _ o O1+0.31 



<^|^^" = 2.21+j;:-^^(stat) ±0.10 (syst) 
</'NC° = 5.21 ± 0.27 (stat) ± 0.38 (syst) . 

An analysis with the constraint of an undistorted ®B energy spectrum 
yielded results which are consistent with previous results in Sec 2. 



Table 2. Systematic uncertainties on fluxes for the spectral shape uncon- 
strained analysis of the salt data set. f denotes CC vs NC anti-correlation. 



Source 


NC uncert. 


CC uncert. 


ES uncert. 




(%) 


{%) 


(%) 


Energy scale 


-3.7,-1-3.6 


-1.0,+1.1 


±1.8 


Energy resolution 


±1.2 


±0.1 


±0.3 


Energy non-linearity 


±0.0 


-0.0,+0.1 


±0.0 


Radial accuracy 


-3.0,+3.5 


-2.6,+2.5 


-2.6,+2.9 


Vertex resolution 


±0.2 


±0.0 


±0.2 


Angular resolution 


±0.2 


±0.2 


±2.4 


Isotropy mean f 


-3.4,+3.1 


-3.4,+2.6 


-0.9,+l.l 


Isotropy resolution 


±0.6 


±0.4 


±0.2 


Radial energy bias 


-2.4,+1.9 


±0.7 


-1.3,+1.2 


Vertex Z accuracy f 


-0.2,+0.3 


±0.1 


±0.1 


Internal background neutrons 


-1.9,+1.8 


±0.0 


±0.0 


Internal background 7's 


±0.1 


±0.1 


±0.0 


Neutron capture 


-2.5, +2. 7 


±0.0 


±0.0 


Clierenkov backgrounds 


-1.1, +0.0 


-1.1, +0.0 


±0.0 


"AV events" 


-0.4,+0.0 


-0.4,+0.0 


±0.0 


Total experimental uncertainty 


-7.3,+7.2 


-4.6,+3.8 


-4.3,+4.5 


Cross section^" 


±1.1 


±1.2 


±0.5 



An analysis to determine the neutrino mixing parameters Am^ and 
tan^ 9 under the assumption of two-flavour MSW neutrino oscillation was 
performed. The ratio fs of the total flirx to the SSM'' value was a free 
parameter, while the total hep flux was fixed at 9.3 x 10^ cm~^ s~^. The 
allowed region of the mixing parameters was determined from a combined 

fit to the D2O and salt data^", and results from other solar neutrino 
experiments^'^'^. The allowed region is shown in Fig. 3(a). If a further 
assumption of CPT invariance is made and the results from the KamLAND 
reactor i7g experiments^ is included, the allowed region is further reduced 
(Fig. 3(b)). The best-fit point with marginalized uncertainties is Am^ = 
7.1^0 3 X 10~^ eV^ and = 32.5tj g degrees, which rejects the hypothesis 
of maximal mixing at a confidence level equivalent to 5.4cr. 
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Figure 2. Distribution of (a) /3i4, (b) cos 00 and (c) volume-weighted radius (p) dis- 
tribution. Also shown are the Monte Carlo predictions for CC, ES, NC + internal and 
external-source neutron events, all scaled to the fit results. The dashed lines represent 
the summed components. All fits were done with To£f>5.5 MeV and Rflt< 550 cm. 
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Figure 3. Global neutrino oscillation contours, (a) Solar global: D2O day and night 
spectra, salt CC, NC, ES fluxes, SK, CI, Ga. The best-fit point is Am^ = 6.5 X 10"^, 
tan^e = 0.40, /a = 1.04, with x^/d.o.f.=70.2/81. (b) Solar global + KamLAND. The 
best-fit point is Am^ = 7.1 X 10"^, tan^ 6» = 0.41, /s = 1.02. In both (a) and (b) the 
®B flux is free and the hep flux is flxcd. 



4. Future Prospects and Summary 

The SNO collaboration is now preparing for the NCD phase of its physics 
program. In this phase, the NC neutrons are detected by the NCD array 
and the CC signals are detected by the PMTs. This separation of the 
detection mechanism will significantly reduce the statistical anti-correlation 
of the signals. In addition, a measurement of the neutrino energy spectrum 
via the CC reaction can be made with higher precision with the reduction 
of NC signal in the Cherenkov light spectrum. 

At the time of this talk, the D2O is being desalinated. All the '^He 
counters have been constructed and are being characterized in the under- 
ground laboratory. There are intense activities in preparing for the NCD 
deployment, and integration of the NCD electronics and data acquisition 
system. 

The SNO experiment has demonstrated conclusively in its pure D2O 
phase that solar p^s transform their flavour while in transit to the Earth. 
A high precision measurement of the total active solar neutrino flux without 
any assumption of the energy dependence of the neutrino flavour transfor- 
mation probability has been made in the salt phase. In the next few months, 
the SNO physics program will enter its NCD phase. It is anticipated that 
the data from this NCD phase, along with data from the previous phases, 
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will significantly enhance our understanding of the properties of neutrinos. 
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